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We examine quantum detection via a Michelson interferometer embedded in a gas with Kerr 
nonlinearity. This nonlinear interferometer is illuminated by pulses of classical light. This strategy 
combines the robustness against practical imperfections of classical light with the improvement 
provided by nonlinear processes. Regarding ultimate quantum limits, we stress that, as a difference 
with linear schemes, the nonlinearity introduces pulse duration as a new variable into play along 
with the energy resources. 
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I. INTRODUCTION 

Precise measurements are crucial in physics since they 
constitute the link between the theory and nature. Ac¬ 
curate measurements can promote or reject a theory. Be¬ 
sides, precise detection and monitoring are fundamental 
for technology and other applications of science. 

A critical contribution of the quantum theory to 
metrology is that quantum fluctuations would limit the 
resolution to some ultimate limits depending on the 
energy resources employed in the process [I|, usually 
counted as the number of particles. 

Since standard metrology is based on linear processes, 
previously known quantum limits heavily depend on an 
implicit assumption of linearity. Thus, a new frontier 
arises if we consider that the signal may be detected 
via nonlinear processes. The key point is that nonlinear 
schemes allow us to reach larger resolution than linear 
ones for the same resources. Moreover, the improvement 
holds even when using probes in classical states. This is 
of much relevance concerning robustness against practi¬ 
cal imperfections, which can be deadly for schemes based 
on nonclassical probe states Quantum nonlinear 

metrology has been studied and proven experimentally 
in very different physical contexts im. In particular, 
this is the case for light propagation in Kerr-type nonlin¬ 
ear media, that has already demonstrated its usefulness 
in the context of precise detection [l3 - [l^ . 

In this paper we present a new feature of quantum 
detection involving nonlinear processes. This is that res¬ 
olution depends not only on the number of photons but 
also on the duration of the pulse (this is both on the 
number of particles and on the rate at which the are em¬ 
ployed) . This is in sharp contrast to linear schemes where 
the duration of the pulse plays no role. In this way a 
new variable appears which may be advantageously used 
to improve detection performance beyond previously ac¬ 
cepted limits. 
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FIG. 1: Scheme of a Michelson interferometer embedded in a 
nonlinear medium. 

II. SCHEME 

To develop this point let us consider signals encoded 
as length variations that can be typically detected via a 
Michelson interferometer schematized in Fig. 1. In the 
absence of a signal the two arms are equal, £2 = ^1 = ^o- 
For simplicity we will consider the signal manifests as 
an anti-correlated length change as = io — x/2, £2 = 
.^0 + ic/2 as it is expected to be the case by the pass of a 
gravitational wave m. 

In order to involve nonlinear effects, we assume that 
the interferometer is embedded in a gas displaying a Kerr 
nonlinearity. The light is made of classical-light pulses of 
mean frequency w, duration r, and cross section A. The 
propagation in the nonlinear medium can be conveniently 
expressed in terms of an intensity-dependent index of re¬ 
fraction, 

n = no + nI = no{l + X^) 1 (1) 

here, n is the nonlinear coefficient, ng is the linear index, 
I is the light intensity, N is the number of photons of 
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each pulse, and x expresses the nonlinear phase shift per 
photon, 
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where this is a definition of x whereas the I versus N 
relation is an estimation good enough for our purposes 
as far as the exact equivalence would require the specific 
spectral distribution of the pulse. 

After Eq. ([Ij, the light propagation within the in¬ 
terferometer is described in the quantum domain by the 
unitary operator U = C/it /2 with 

= g,=N, + ^NI (3) 


where rj is the quantum efficiency of the detectors and a 
is the variance of the random relative phase. 

In order to be detected, the signal-induced shift {M) 
must be larger than the background quantum noise AM 
at X = 0. Taking into account the noise sources com¬ 
mented on above we get 

(AMf = riN + rfN'^a^ + r^NNt, (5) 

where Nt is the mean number of the undesired thermal 
photons. Since the nonlinearity complicates the noise 
analysis we postpone the details of the calculation to an 
Appendix to not interrupt the analysis. 

We can estimate the resolution of the detection of x 
via the noise-to-signal ratio as 


with ipj = kij where k = uqujIc is the wave number, 
and Nj = a^aj represents the photon-number operator 
in each arm j = 1,2. We further assume that the signal 
induces an arm length difference x small compared with 
the length of the pulse cr. Any other changes produced 
in the optical constants of the media are assumed to lead 
to negligible effects. 

We illuminate the interferometer just by one of the 
input ports (as usual the vacuum is at the other in¬ 
put) by a classical-like pure coherent state |a) with a 
mean number of photons |ap = N ^ 1 [exactly the 
same results are obtained if the probe state before the 
input beam splitter is in the phase-averaged mixed state 
d9\ae^^){ae^^\/(2 -k)]. The light state in the internal 
modes ai _2 of the interferometer leaving the 50% beam 
splitter can be expressed as the product of coherent states 
|a/-\/ 2 )i|a/-\/ 2 ) 2 , whereas the light state reaching the 
beam splitter after propagation within the interferometer 
is UiU 2 \a/V 2 )i\a/^/ 2 ) 2 . 

The measurement is carried out by registering the dif¬ 
ference in the number of photons recorded by two de¬ 
tectors at the output ports of the interferometer. The 
corresponding operator can be expressed in terms of 
the internal modes of the interferometer 01^2 as M = 
i{alai — a|a2). 

The sources of technical noise will be taken into ac¬ 
count by their most typical consequences such as phase 
randomization (e.g., caused by fluctuations in the optical 
properties of the medium), thermalization (e.g., thermal 
photons coming from undesired residual sources), as well 
as the usual finite quantum efficiency of the detectors. 


III. SIGNAL DETECTION AND 
UNCERTAINTY 
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leading to 


^2 1 + riNa'^ + Nt 

(Ax) = - 7T. 

r]k^N{l + xN/2)^ 


(7) 


This holds provided that yA/ca: <C 1, iV, and 

CT <C 1. Moreover, we have assumed xfc£o — 27rm for 
integer m. Otherwise, without this last condition the 
nonlinearity would deeply disturb the probe state by pro¬ 
ducing coherent superpositions of distinguishable states 
[l8l| . and the standard interferometric measurement M 
would become useless requiring more advanced detection 
strategies beyond the scope of the present analysis (see 
the Appendix for further details). The nonlinearity will 
have a noticeable effect for yA ^ 1, which is compatible 
with the above assumption yAfca: <C 1 provided that kx 
is small enough kx <^1. 

This might be compared with the case when the non¬ 
linear medium is absent y = 0 
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and we have further assumed that in such a case prop¬ 
agation occurs in vacuum and the phase randomization 
can be safely neglected cr = 0. 

In the ideal case that the phase randomization and 
thermal effects might be ignored a = Nt = 0 and yA 1 
we get the following improvement of the nonlinear versus 
the linear scheme. 


The signal x produces a phase shift that alters the 
statistics of the observed M, shifting its mean value (see 
the Appendix for details). 
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We recall that the duration of the pulses is embedded in 
the nonlinear phase shift per photon y in Eq. and so, 
the lesser r, the larger y, and the larger the resolution. 









3 


IV. DISCUSSION V. CONCLUSIONS 


We can roughly estimate the amount of noise reduction 
with parameters within the reach of current technology. 
For the sake of simplicity and to fix the main ideas let us 
first consider the ideal case where the effect of technical 
noise is negligible rjNa^ + TVj <C 1 and ry ~ 1. Regarding 
numerical values we can address two extreme situations: 
standard natural nonlinearities, and giant nonlinearities 
achieved via atomic coherence. Throughout we will as¬ 
sume that the index in darkness is on the order of unity 
no ~ 1 . 

The typical natural nonlinearities in gases can be on 
the order of h ~ 10“^^cm^/W [l^. As to the pulse 
parameters let us assume a pulse duration of r ~ 1 ps, 
light power of P ~ 1 PW, and beam cross section of A ~ 
10“®m^, which leads in the visible spectrum to V ~ 10^^ 
photons per pulse and a nonlinear shift per photon of 
X ~ 10“^®, so that 

Ax ~ 10“® ~ 10“^^m. (10) 

The condition x^kx <C 1 means the following condition 
on the signal x 10“^°m. On the other hand, the condi¬ 
tion xkio — 27rm leads to an extremely large interferome¬ 
ter even for m = 1 since in such a case £o ~ lO^^m. Thus 
the m = 0 situation should be addressed as suggested in 
the Appendix. 

Things are completely different if we consider the gi¬ 
ant nonlinearities achieved via electromagnetically in¬ 
duced transparency, leading to Kerr coefficients on the 
order of h ~ 10“^cm^/W as reported in Ref. [2^ and 
similarly large values in other configurations, such as 
h ~ 10“®cm^/W in Ref. [2l|. Such large values al¬ 
lows for alleviating the requirements on the light probe 
state. For example we may have t ~ 100 ps, P ~ 1 MW, 
and A ~ 10“®m^, which leads in the visible spectrum to 
N ~ 10^'* and x — 10“® so that 

Ax ~ 10“® ~ 10“^°m. (11) 

The condition xk£o — 27rm leads to a more practical 
interferometer with £o — 100 m, whereas the condition 
xNkx <C 1 implies x <C 10“^®m. Thus a detectable 
signal should be composed in the range 10“^®m ^ x ^ 
10“^°m. This fits perfectly well with the expected signals 
due to the pass of a gravitational wave in a 100-m-long 
interferometer, which are 10“^®m ^ x ^ 10“^°m [^ . 
Notably, smaller r and/or A, such as the beam-size val¬ 
ues reached in Ref. |23l| with current technology, may 
lead even to room-size interferometers with similar per¬ 
formance. 

Finally, we may estimate the maximum effect of imper¬ 
fections so that the good effects of nonlinearity are not 
spoiled. The condition we are looking for is derived from 
rjNcr^ + Nt x^k^^ that is satisfied if, roughly speaking, 
(7 <C X\/^/'n ^ For the natural nonlin¬ 

earity and ry ~ 1 we get cr <C 10“®, Nt 10®, whereas 
for giant nonlinearities we get much less limiting bounds, 
a < 10 - 1 , Nt < 1012 . 


Summarizing, nonlinearity not only can improve res¬ 
olution beyond linear limits, but also introduces a new 
variable into play. The signal uncertainty depends not 
only on the number of probe photons V, but also on the 
duration of the pulse r through the nonlinear effect per 
photon X in Eq. ©• This is because optical nonlinearity 
is sensible to light intensity rather than just energy or 
photon number. In particular, after Eq. 0 we may con¬ 
jecture an optimum ultimate quantum limit (that would 
require nonclassical probes to be reached) scaling as 

tA\^ 

( 12 ) 

in terms of the probe free parameters, where A is the 
wavelength. This result may be particularly useful for ex¬ 
ample in situations of frequent monitoring where small 
pulse durations and large repetition rate of the inter¬ 
rogating pulse may be of interest. In this regard, the 
availability to obtain large beam intensities by shorten¬ 
ing pulses seems a more feasible condition than increas¬ 
ing energy resources as required in usual linear quantum 
metrology. 
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Appendix A: Calculus 

After the relation aF{a^a) = F(ala-|-l)a valid for any 
F we get that 

(Al) 

where Lpj = k£j, Zj = (pjxl‘2, and U = U 1 U 2 is the 
global transformation. When evaluating the mean values 
of a\a 2 and its Hermitian conjugate on coherent states 
1/3) with /3 = aj^, we will get expressions of the form 

(/3|C/1'a^-C/|/3) = e*‘^^'e*^^(/3|e*2^^-^^a^-|/3), (A2) 

that can be easily evaluated taking into account that 

(/3|e*2^^'^^aj|/3) = ^(/3|/3e*2^0 = (A3) 

Besides the finite quantum efficiency of the detectors 
we will consider some further typical forms of practical 
noise, such as thermalization and phase randomization. 
These common noise forms can have different physical 
origins, such as fluctuations of the optical properties of 
the medium, random variations of the complex amplitude 
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from pulse to pulse, and so on. They can be addressed 
at once by performing the replacements, 

a\a 2 -t + \/l - (v^a2 + v^l - 7762) , 

(A4) 

where ^ is a random phase that we will assume to be 
Gaussian distributed with zero mean and variance cr^, 
r] is the quantum efficiency in the detection, and bj^s 
are uncorrelated field modes in thermal states with (1 — 


r]){b\bi) = (1 — f]){blb 2 ) = Nt/2 with Nt N, and 
(bi) = ( 62 ) = {b\b 2 ) = 0 . 

1. Mean value 

Taking all this into account the mean value of M can 
be obtained after a long but straightforward calculation 
as 





Before the (/) integration several natural considerations 
seem in order to get simpler and meaningful expressions. 
A required condition is that the factor in the real ex¬ 
ponential should be close to zero, otherwise the final 
uncertainty Aa; would increase exponentially with N. 
This is because the uncertainty AM will contain always 
a photon-counting noise term independent of the arm 
lengths. Thus we have to consider that in the absence 
of signal zq = tott, where m is any integer. This may be 
achieved by properly adjusting the fixed arm length 
depending on x- Alternatively we may consider that the 
Kerr transformation induced by the fixed length may 
be compensated by another Kerr transformation with a 
nonlinear susceptibility of the opposite sign in propaga¬ 
tion conditions insensitive to the signal value x. An alter¬ 
native approach may follow the strategy in Ref. [ 2 ^ by 
comparing outputs for two consecutive pulses experienc¬ 
ing alternatively linear and nonlinear transformations. 

Thus, considering that the signal induces a very small 
variation in Zj around zo — mw we have 


(M) ~ sin 


(/> + kx 



(A6) 


which can be obtained after Eq. (IA5I) by a series expan¬ 
sion of the harmonic functions within the exponential 
and the sine function where we have also neglected the 
Z 2 — zi term not multiplied by N. Carrying out the (j) 
integration over a Gaussian distribution with zero mean 
and variance we get 


kx{l + XY 

(A7) 

Finally, considering typical values for the variables in¬ 
volved the real exponentials can be safely approximated 
by unity. Moreover, the expected signals are small 
enough x^kx <C 1 so that there is a linear relationship 
between M and x, which is usually an implicit assump- 


(M) ~ ijNe 




sm 



N j 

- <pi + Z 2 - zi + Y [sin ( 2 ^ 2 ) - sin (2zi)] ^ . (A5) 


tion leading to Eq. m, 


(M) ~ r]N 


kx 



(A8) 


This is compatible with the fact that yA can be very 
large having observable effects. This is because it deter¬ 
mines the value of the x derivative in the denominator 
of Eq. ([5]) and such a derivative need not be small. On 
the contrary, the best situation holds when xA is large 
enough to imply a noticeable reduction in signal uncer¬ 
tainty as analyzed in detail in Secs. Ill and IV. 


2. Uncertainty 


Next we address the evaluation of (AM)^ at ^2 = ^1 = 
£ 0 , this is x = 0 so that (M) = 0 and (AM)^ = (M^). 
After Eq. (jA4l) the effect of thermalization and finite 
efficiency means that M^ should be replaced by 


- rf) 


+ l) Ai + {2b\bi + l) N 2 


(A9) 

where to avoid confusion we denote by Mq when evaluat¬ 
ing M in the noiseless case. Other terms lead to null or 
negligible contributions. It is worth noting that the last 
term is not affected by the Uj transformation nor by the 
random phase. 

Then we can compute (Mq), where 


= 2A1A2 +N1+N2- of ai - al^al (AlO) 


The first terms depending just on Nj are invariant under 
the transformations Uj whereas for the remaining two 
terms we can use the fact that a^F(a^a) = F{a^a + 2)a? 
and then proceed as above, to get, before phase random¬ 
ization. 


o A^ A^ 

{Ml) = — + N-—cos{2ct>\ 


(All) 
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and after the random-phase average, 

/V2 Ar2 2 

(Mp^) = —-H7V-—~7V-tcr2Ar2^ (A12) 

where the approximation holds for a I. Finally, col¬ 
lecting the contributions in Eqs. (IA9I) and (IA12I) we fi¬ 


nally get 

(AM)^ ~ rjN + + rjNNt, 

leading to Eq. ([7]). 


(A13) 
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